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Abstract 

The production rates of D** mesons in charm and bottom events at centre-of-mass energies 
of about 91 GeV and the partial width of primary cc pairs in hadronic decays have 
been measured at LEP using almost 4.4 million hadronic decays collected with the 
OPAL detector between 1990 and 1995. Using a combination of several charm quark 
tagging methods based on fully and partially reconstructed D** mesons, and a bottom 
tag based on identified muons and electrons, the hadronisation fractions of charm and 
bottom quarks into D** mesons have been found to be 

f(b ^ D*+X) = 0.173 ±0.016 ±0.012 and f (c ^ D*+X) = 0.222 ± 0.014 ± 0.014 . 

The fraction of cc events in hadronic Z^ decays, Pcc/Phad = T{Z^ cc)/T{Z^ hadrons), 
is determined to be 

Tcc/rhad = 0.180 ± 0.011 ± 0.012 ± 0.006 . 

In all cases the first error is statistical, and the second one systematic. The last error 
quoted for Pcc/Phad is due to external branching ratios. 
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1 Introduction 



The production of heavy quarks in the decay of the Z° boson and their hadronisation have been 
the subject of considerable interest over the last few years. In particular the fractions with which 
the Z° boson decays into quark pairs of flavour q have been studied extensively in Z° bb 
decays in Z° — > cc decays and in light flavour events 0. The fraction of bb events 

in Z° decays has been measured with very good precision. To achieve this goal, very efficient 
and pure bottom tagging methods have been developed, resulting in samples of events that are 
nearly free of non-bottom backgrounds. Significantly fewer and less precise measurements exist 
of the equivalent quantity for cc events or for light fiavour events. In particular the selection of 
a pure cc sample has met with many difficulties, and the efficiencies and purities achieved by 
charm tags are inferior to those for bottom tags. The reason for this is that charmed hadrons 
are lighter and shorter lived than bottom hadrons, and are similar enough to most light hadrons 
to make a separation very difficult. However, the precise knowledge of the partial widths for 
different flavours constitutes an important test of the predictions of the Standard Model, since 
in lowest-order Born approximation the partial Z° decay width to qq, Fqq, is related to the 
coupling constants of the vector and axial vector current, and g'j^: 

^m = N^^mf + {9lf)- (1) 

Here is the Fermi decay constant and mz the TP mass. The factor = 3 denotes the 
number of colours. Higher order electroweak and QCD corrections to the propagator and qq 
vertex that modify Fqq essentially cancel in the ratio Fqq/Fhad except in the case of Z" — >• bb, 
where a small dependence on the Higgs mass and on the precise value of the top quark mass 
is introduced. The ratio Fqq/Fhad therefore is the preferred measurable quantity, for which 
precise predictions exist in the context of the Standard Model for the quark flavours u,d,s and 
c, almost independent of unknown quantities. 

In this paper a measurement of the fraction of primary cc pairs produced in the decays 
of bosons is presented. At the same time the hadronisation fractions f (c — D*+X) and 
f (b — D*'''X) are measured. The analysis is based on the identiflcation of charged 0*=*= mesons, 
electrons and muons. 

The hadronisation fractions f (c ^ D*+X) and f (b D*+X) are measured using a double 
tagging technique. To determine f (c ^ D*+X), charged 0*=*= mesons are sought in both event 
hemispheres]]. The hadronisation fraction f (b ^ D*+X) is determined in events tagged by a 
hard lepton in one hemisphere and a D** in the other hemisphere. Comparing the number of 
such double tagged events with the number of singly reconstructed 0*=*= mesons or leptons, the 
hadronisation fractions can be extracted with minimal model dependence and without explicit 
knowledge of the D** or lepton reconstruction efficiencies. 

The ratio of the charm partial width to the total hadronic width, Fcc/Fhad, is determined 
from the hadronisation fraction f (c — D*+X) and from a measurement of the total production 
rate of D** mesons in Z° cc events, Fcc/Fhad'f (c — > D*+X). This rate is measured in this paper 
using a particularly well understood D*''= decay mode, the decayf] D*+ — > D^vr"^, D° — > K'tt"*". 

Both the measurement of the hadronisation fraction and the measurement of Fcc/Fhad rely 
heavily on the reconstruction of D*+ mesons using two different techniques. Therefore the 

^The plane separating the two hemispheres in an event is defined perpendicular to the thrust axis of the 
event. 

^ Charge conjugation is assumed throughout this paper. 
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discussions in the first part of the paper concentrate on these technical aspects of the analysis. 
In the first technique described in section | the D*+ mesons are identified in a number of different 
decay channels by reconstructing all charged decay products. Since a significant contribution 
to the D*+ sample is from bottom hadron decays, a method has been developed to separate 
the different sources and is discussed in some detail. The second method of D*+ reconstruction 
is described next. It is a much more inclusive method, where only the pion in the decay 
D*+ —>■ D^TT~^ is used as the tag for the D*+. In the last part of the section the tagging of 
Z" — > bb events using leptons is summarised. 

In the second part of the paper the different measurements are presented. In section ^ the 
determination of the total rate of 0*+ production in hadronic Z° cc decays is described. This 
is followed in section |] by a presentation of the double tagging technique used to measure the 
hadronisation fractions for both bottom and charm events. Finally the results are combined 
to derive the relative partial width Fcc/Fhad- The results reported in this paper supersede the 
ones given in and complement the measurement of Fcc/Fhad reported in 0. 



2 Analysis Principle 

The main goal of this analysis is the measurement of the hadronisation fractions f (c — » D*'''X) 
and f (b — s> D*+X) and of Fcc/Fhad- A double tagging technique is used to minimise model 
dependencies. However, because charm tags are rather inefficient a full double tag, where the 
same tag is applied to both hemispheres of the event, cannot be used. Instead two different 
charm tags are applied, one, which is pure, but has a comparatively small efficiency; and the 
other, which is more efficient, but less pure. The general strategy for the measurement of 
the hadronisation fractions is that a charm or bottom enriched sample is selected by applying 
the high purity charm or bottom tag to one hemisphere of the event, and then searching for 
D*+ mesons in the opposite hemisphere using the more efficient, less pure tag. Neglecting for 
simplicity any background from other fiavours, the number of events of fiavour q is given by 

iVtagi ~ ^ f (q - D*+X) etagi , (2) 

J- had 

where etagi is the efficiency to select an event of fiavour q using the pure tag. The number of 
events, where a D*"*" mesons is simultaneously identified in the second hemisphere, is given by 

iVtagltag2 ~ -P^f(q^D*+X)etaglXf(q^D*+X)etag2 
-Lhad 

= iVtagi xf(q-.D*+X)etag2 . (3) 

Here etag2 is the efficiency to tag a D*+ mesons using the second, efficient, tag in the fiavour 
tagged sample. From the ratio of the number of double tagged events to the number of single 
tagged events the hadronisation fraction f (q — > D*+X) can be measured essentially without 
further assumptions or inputs. In this analysis the high purity fiavour tags are an exclusive 
D*+ tag for Z° — > cc events, and a lepton tag for Z° — >• bb events. The D*+ tag, applied to the 
sample of fiavour tagged events, is based on a very inclusive method of identifying D*+ mesons 
using only the pion from the decay D*+ D°7r+. 

Significant backgrounds however exist from other than the desired fiavours. In addition the 
efficiency to find a D*+ meson in the fiavour tagged sample is not independent of the fiavour 
tag itself. Background is particularly important in the case where the fiavour tag is the charm 
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tag. A significant part of tlie sample of D*"*" mesons originates from TP bb events, and 
also non-negligible contributions from combinatorial background events are found. The charm 
tagged sample is selected by fully reconstructing D*+ mesons that decay into a particular final 
state K~Y with a branching ratio E = 6(0*+ — »• D^tt"*") B(D'' — > K~Y). The number of such 
events, iVD*+ , is given by 

iVD.+ = 2iVhad ■ E f (q ^ D*+X) B e^,+ ) + iV^^^ (4) 

q=b,c -"-had 

Here, A^had is the number of hadronic Z° decays used, Fq^/Fhad is the relative partial width for a 
Z° to decay into a quark-antiquark pair of flavour q, /q is the fraction of events with flavour q 
in the sample, 6^*+ is the efficiency to reconstruct a D*"*" meson in a q — D*+ hemisphere, and 
iVbgd is the number of background events in the sample. The fractions /q satisfy the condition 

/c + /b = 1. 

In this sample of flavour tagged events D*+ mesons are sought in the opposite hemisphere 
using the inclusive D*+ reconstruction. Background in the sample is reduced by requiring that 
the two D*"*" candidates have opposite charges. The contribution from b —>■ D*+ decays is further 
reduced since some of the bottom hadrons will have mixed before decaying into a D*+ meson. 
The number of double tagged events is therefore given by 

where i3* = B(D*+ — » D^tt^), e^*^ is the efficiency for tagging a D*+ meson using the inclusive 
tagging method in a Z° — >■ qq event, when a D* meson has been identified in the opposite 
hemisphere, and Xefr is an effective B°— B mixing parameter applicable to the selected sample 
of events. It is interesting to note that eq. |^ does not depend on the efficiency of the high purity 
flavour tag, but only on the efficiency of reconstructing a D*+ meson inclusively. 

The number of double tagged events given in equation]^ is proportional to both f (c ^ D*+X) 
and to f (b — > D*+X). To separate the components the analysis is done twice, once as shown 
in eq. |^ in a charm enriched sample, tagged by the presence of D*+ mesons, and once in a 
bottom enriched sample, selected through hard leptons. The latter analysis is mostly sensitive 
to f (b — * D*'*'X), the former to f (c ^ D*+X). By fitting the two samples simultaneously both 
hadronisation fractions are determined. 

The charm partial width, Fcc/Fhad, is determined from the measured hadronisation fraction 
f (c — i> D*+X) and from the total rate with which D*"*" mesons are produced in Z° — cc decays, 
Fcc/Fhad f (c ^ D*+X). This is measured using a particularly well understood D*+ decay mode, 
D*+ —>■ D°7r^,D° K^TT^, which facilitates the source separation and the background sub- 
traction. The only quantity which is not measured in this analysis, but has to be taken from 
external sources, is the branching ratio B(D*+ D^tt^). 

3 The OPAL Detector and Event Selection 

A detailed description of the OPAL detector can be found elsewhere []T0[. This analysis relies 
heavily on the precise reconstruction of charged particle tracks and primary and secondary 
vertices in the event. This is achieved using a combination of two layers of a high precision 
silicon micro- vertex detector, installed nearest to the primary interaction point, and a system 
of large-volume gas-filled drift chambers which combine excellent spatial resolution with very 
good particle identification capabilities via the measurement of the specific energy loss of tracks. 



/, f (c D*+X) e^.^ + /b (1-Xeff) f (b - D*+X) 6^,^ -B^+N^t , (5) 
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The whole central tracking system is immersed in a magnetic field of 0.435 T, oriented along 
the direction of the electron beam. These central tracking detectors are surrounded by both 
electromagnetic and hadronic calorimeters with good energy resolution, providing nearly her- 
metic coverage over the full solid angle, and by a system of muon chambers on the outside of 
the detector. 

Hadronic Z° decays are selected based on the number of reconstructed charged tracks and the 
energy deposited in the calorimeter []TT|. The total hadronic event selection efficiency is found 
to be (98.7 ±0.4)%. The selection slightly changes the fiavour composition of the sample. This 
flavour bias is found to be less than 0.1% |Q. The analysis uses an initial sample of 4374410 
hadronic decays of the Z° collected with the OPAL detector between 1990 and 1995. 

Jets are reconstructed in the events by the cone jet finder |]12| with a cone radius, -Rcone, set 
to 0.7, and a minimum cone energy of at least 5 GeV. Events are only accepted if at least two 
jets are reconstructed. To ensure that the event is mostly contained in the sensitive detector 
volume, the absolute value of the cosine of the polar angle of the thrust axis with respect to 
the beam direction, | cos 6'thrust | , has to be smaller than 0.9. 

Tracks are used in the reconstruction if they pass loose track quality cuts requiring \do\ < 0.5 
cm, \zo\ < 20 cm, p^y > 0.250 GeV and ^cj > 40. Here do is the distance of closest approach 
between the primary vertex and the track measured in the plane perpendicular to the beam, 
zo the distance along the beam at this point, p^y the momentum in the plane perpendicular 
to the beam, and ncj the number of hits on the track which are reconstructed in the main 
tracking chamber. The primary vertex in a collision is reconstructed from the charged tracks 
in the event and constrained by the known average position and spread of the e"'"e~ interaction 
point. 

Hadronic decays of the Z° have been simulated using the JETSET 7.4 Monte Carlo model 



||13|| with parameters tuned to the data |T^. A sample of approximately 10 million simulated 
events was available for this analysis. In all samples heavy quark fragmentation was imple- 
mented using the Peterson model with fragmentation parameters determined from LEP 
data [16]. All samples have been passed through a detailed simulation of the OPAL detector [l7 
before being analysed using the same programs as for data. 



4 Heavy Flavour Tagging Techniques 

Three different tagging techniques are employed to identify Z° cc and Z° — > bb events. The 
charm tags are based on the exclusive reconstruction of charged D*""" mesons (called "exclusive 
tag" in the following) in five different decay chains, or on a more inclusive D*+ reconstruction 
(called "inclusive tag" in the following). Bottom events are identified through the presence 
of an electron or a muon with large momentum and large transverse momentum relative to 
the direction of the jet containing the lepton. In this section the different tagging methods 
are described in some detail. Particular emphasis is placed on the method used to separate 
the contributions to the samples tagged by the different methods, where large backgrounds 
are present, and on the systematic errors connected with this source separation method. Very 
similar techniques have been used in previous OPAL publications for charm tags 0,|l5] and the 
bottom tag using leptons @,|19 . 
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4.1 The Exclusive Charm Tag 
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The exclusive charm tag is based on the reconstruction of charged D*+ mesons in five different 
decay channels: 

D*+ D^TT + 

U K-TT+ "3-prong" 

' — > K~e~^h'e "electron' 

' — > ¥<r ji^u^ "muoh 

U K-tt+ttO "satellite" 

' — > K^7r+7r^7r+ "5-prong" 

In the following the electron and the muon channel are collectively referred to as "semileptonic" . 
No attempt is made to reconstruct the 7r° in the satellite channel, or the neutrino in the two 
semileptonic channels. Electrons are identified based on their energy loss in the jet chamber and 
the energy deposition in the electromagnetic calorimeter. An artificial neural network trained 
on simulated events is used to perform the selection Muon candidates are identified by 
associating tracks found in the central tracking system with tracks in the outer muon chambers 

In each channel a candidate is formed by combining an appropriate number of tracks, 
corresponding to the number of charged decay products in the chosen decay mode, assigning 
one to be a kaon, the rest to be pions or leptons, and calculating the invariant mass, Mq, of 
the set of tracks. Candidates are selected if the tracks assigned to the decay products have the 
correct charges, and if the reconstructed mass lies within the expected range, defined by the 
mass resolution in the different channels. The exact values are given in table |l|. After adding a 
further track as a candidate for the pion in the D*+ decay the combined mass, M^,, is calculated 
and the candidate is selected if the mass difference AM = M^, — Mq is within given limits. Note 
that in cases where not all particles from the particular decay are reconstructed, the masses Mq 
or do not correspond to the physical particle masses of the D° or D*+ mesons, respectively. 

For candidates with xd*+ = E^^+ /E\,ea.m < 0.5 the particle identification power of the OPAL 
detector is used to enrich the sample in true kaons. A probability Wf^j^^, that the difference 
between the measured specific energy loss, dE/dx, determined for a track of momentum p, and 
the d-E/dx value expected at that momentum for a kaon, is compatible with the kaon particle 
hypothesis, is calculated. A candidate track has to have a probability of at least 2% to be 
accepted as a kaon candidate. To ensure a reliable dE/dx measurement the number of charge 
measurements used in the dE/dx calculation, n^E/Axi has to be at least 20. 

Background in the sample is further reduced by cutting on the helicity angle 9* , measured 
between the direction of the D° candidate in the laboratory frame and the direction of the 
kaon in the rest frame of the D° candidate. True kaons from D" decays are expected to be 
isotropically distributed in cos 6'*, while background displays pronounced peaks at cos 6'* = — 1 
and, particularly at low xd*+, at cos^^* = +1. 

To avoid multiple counting of events if more than one D*+ candidate is found, only one 
candidate per event is accepted. If more than one candidate is found per event a hierarchy is 
used to select the best one according to the signal purity. A 3-prong decay is preferred over 
a semileptonic one, which in turn is preferred over a satellite, and a 5-prong is selected last. 
For the semileptonic channel an electron is preferred over a muon candidate. If more than 
one candidate is found within one channel, the one with Mq closest to its nominal value of 
Mdo = 1.864 GeV ^ (1.6 GeV for the satellite) is selected. 

A detailed list of the cuts used in the different channels is given in table 0. The reconstructed 
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Jb 1 CLlitlC 


3-prong 


semileptonic 


satellite 


U JJi Clig, 






0.4-1.0 


0.4-1.0 


0.4-1.0 


0.5-1.0 


Mq [ GeVJ 


mil 


1.79-1.94 


1.20-1.80 


1.41-1.77 


1. (9-1.94 


AM [ GeV] 


full 


0.142-0.149 


0.140-0.162 


0.141-0.151 


0.142-0.149 




< 0.5 


> 0.02 




^dE/dx 


< 0.5 


20 




COS 9* 


< 0.5 
> 0.5 




-0.8-0.8 
-0.9-1.0 




-0.9-1.0 



Table 1: List of selection cuts used in the D*"*" reconstruction. Note that both the scaled energy 
x-£)*+ and the mass Mq are effective quantities, calculated from the reconstructed tracks only. The 
exact meaning of the different quantities is explained in the text. 

invariant mass spectra in all channels exhibit characteristic peaks at AM = 0.145 GeV, close 
to the kinematic threshold of AM = m^r. Clear signals are visible in all five channels, as shown 
in figure 

4.1.1 Partially Reconstructed D*+ Mesons 

A significant fraction of the sample of selected D*+ mesons are only partially reconstructed. 
This is particularly obvious for events where neutral decay products are not identified, as is the 
case for the satellite or the semileptonic channels. Even in decays where all decay products are 
charged particles a fraction of events is present in the sample which are real D*+ decays, where 
however one or more tracks have been wrongly identified, or are missed completely. Such events 
are called "partially reconstructed D*"*" mesons" if the slow pion of the D*+ — > D^vr^ decay has 
been correctly found. These partially reconstructed D*""" mesons produce an enhancement in 
the AM spectrum very similar to the true signal. Only very few of these events are present 
in the 3-prong sample. They are much more important in the 5-prong tagged events, where a 
clear tail is visible in the AM distribution for values above 0.145 GeV (see figure |I|(d)). Since 
such events originate from D*+ decays, they can still be used in the flavour tagging part of the 
analysis. 

4.1.2 Combinatorial Background Estimation 

The dominant background source is random combinations of tracks that pass the applied cuts. 
Only this combinatorial background component is considered background for the flavour tag- 
ging, and a method has been developed to subtract only this component from the sample of 
tagged events. The combinatorial background component is described by an estimator con- 
structed entirely from data, optimised to exclude partially reconstructed D*+ candidates. It is 
constructed using a hemisphere mixing technique first introduced in [^. The candidate for 
the pion in the D*""" D°7r+ is taken from the opposite hemisphere relative to the rest of the 
candidate, and reflected through the origin, before being used in the calculation of the invariant 
mass. No requirements are placed on the charge of the D*""" background candidate, except that 
the total charge should be ±1. The resulting distribution is used to define the shape of the 
background in AM. This method ensures that no true pions from the D*+ D^Tr^ decay are 
included in the estimator, and that the background shape does not exhibit any peak in the 
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Figure 1: Distributions of the difference AM = M^, — Mq reconstructed in the four different D*"*" 
channels. The arrows indicate the selected signal region, (a) 3-prong decay mode, (b) the two semilep- 
tonic modes combined, (c) the satellite decay mode, and (d) the 5-prong decay mode. The points with 
error bars are the signal candidates. Superimposed in each case (line histogram) are the background 
estimator distributions, normalised to the upper sidebands in AM. 
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decay channel 


xy)*+ range 


-^cand 


Nhgd 


3-prong 


0.4-1.0 


4649 


1034±28 


semileptonic 


0.4-1.0 


2485 


587±23 


satellite 


0.4-1.0 


10086 


4537±64 


5-prong 


0.5 - 1.0 


9785 


5208±64 


total 




27005 


11366±107 



Table 2: Number of observed candidates in the signal region and the estimated number of these 
which are background events. The error quoted for the background is the statistical error of the 
background sample, and does not contain systematic effects. 



interesting AM region. The background distribution thus obtained is normalised to the candi- 
date AM distribution in the range 0.18 GeV < AM < 0.20 GeV (0.19 GeV < AM < 0.22 GeV 
in the semileptonic channels). Monte Carlo studies have shown that this "reflected pion" esti- 
mator reliably models the shape of the combinatorial background in the sample. The number 
of candidates and the estimated number of background events are given for all channels in 
table g. 

4.1.3 Flavour Composition of the Exclusive 0*=*= Sample 

The D*+ mesons contained in the tagged sample originate mostly in charm and bottom events, 
with a small contribution from events where a gluon splits into a pair of charm quarks. The 
latter is highly suppressed because of the high xd*+ cut applied to the sample of selected events. 
Even more suppressed is the production of D*+ mesons from gluon splitting events into pairs of 
bottom quarks because of the large mass of the bottom quark. In this section the determination 
of the composition of the sample is described. 

The composition of the tagged events is determined by applying three different bottom tags 
to the samples^, and combining the results. Assuming for simplicity that no background from 
light flavours is present in the sample, the number of events tagged by a particular bottom tag 
is given by 

A^b-tag = (/b n + fc Vc) Ar,,,d • (6) 

Here /b = 1 — /c is the bottom (charm) fraction in the sample, and the Vh,c are the probabilities 
that an event in which a D*+ meson has been identifled is also tagged by the bottom tag. If 
these tagging probabilities are known, the bottom fraction can be calculated. The bottom 
tags are applied on a jet basis, both in the jet containing the exclusively reconstructed D*+ 
candidate (called the "D-jet"), and in the remaining highest energetic jet in the event (called 
the "secondary jet"). 

In practise additional backgrounds are present in the sample. The number of background 
events tagged with the different bottom tags and included in the sample is measured from data 
in independent background samples. For this a background sample for each tagging technique 
is prepared, the total number of events in the background sample is normalised to the number 
of background events measured, and the bottom tag is applied to this sample. 

•^The same method has been used in M, where additional details may be found. 
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The three different bottom tagging techniques are based on hfetime information, on jet- 
shapes and on hemisphere charge information. The first two have been used in earher OPAL 
pubhcations [0,|18[, and are only briefly reviewed. The last one will be covered in more detail. 

Lifetime information is reconstructed in both jets used in this analysis. Vertices are recon- 



structed inclusively as in and a decay length significance d/a is calculated, where d is the 
distance between the primary and the secondary vertex, constrained by the jet direction, and 
a its error. Bottom events are identified by their large decay length significance values. The 
shape of the combinatorial background is estimated using the reflected pion technique discussed 
in section [4.1.2|. The background estimator distribution is normalised to the sidebands in AM, 



and is subtracted from the candidate distributions. 

Jet-shape information is used in the jet opposite the D-jet. The shapes are measured by 
a set of seven jet shape variables, which are defined in 0], and are combined using a neural 
net technique into one tag. The combinatorial background is estimated in data using a wrong 
charge technique, where background events are identified by the presence of a candidate with 
an unphysical charge combination of the decay products. 

The third method uses the observation that the charge of the primary quark can be measured 
on a statistical basis using the hemisphere charge. Since the correlations between the charges 
of the D*+ mesons and the sign of the charge of the primary quark are opposite for bottom 
and for charm quarks, measuring the D*+ charge and the primary quark charge in the opposite 
hemisphere provides some separation between bottom and charm events. The hemisphere 
charge is determined from all tracks in the hemisphere according to 

^ _ X^i \Pi\ Qi /r,\ 

EM^ ' ^ ^ 

where i runs over all tracks in a hemisphere, pi is the momentum component along the thrust 
axis of track i in the hemisphere, and qi is its charge. The exponent k is a weighting factor which 
has been optimised using Monte Carlo simulation to be 0.4 for the purpose of flavour separation. 
Similar to the other two methods described a tagging efficiency Pq, q = b, c is determined from 
Monte Carlo. The shape of the background in the hemisphere charge is estimated from events 
tagged in sidebands of the AM distributions, AM > 0.18 GeV, in the D*+ sample. 

In all three cases the tagging probabilities are taken from the Monte Carlo simulation. The 
final fit for the flavour composition is performed simultaneously with the information from all 
three methods. It is done separately for each exclusive channel considered, and in bins of the 
scaled energy xd*+ of the candidate. The most significant contribution to the separation comes 
from the decay length significance analysis, which contributes with a weight of 0.41 from the 
D*+ hemisphere analysis, and 0.27 from the opposite hemisphere. The jet-shape analysis enters 
with a relative weight of 0.21, the hemisphere charge with 0.11. Distributions of the tagging 
variables are shown in figure |^. Combining all exclusive channels the charm fraction in the D*+ 
tagged sample, for the range of detailed in table is determined to be 

= 0.774 ± 0.008 , (8) 

where the error quoted is purely statistical. 

4.1.4 Contribution from g — > cc 

Small contributions to the signal are expected from the splitting of a gluon into a pair of charm 
quarks. This rate has been measured in [^|2^, where the multiplicity of cc pairs produced in 
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Figure 2: Distributions of the different tag- variables used in the flavour separation. Shown are 
the data distribution after background subtraction (points with error bars), the equivalent Monte 
Carlo distribution for all candidates (open histogram) and the predicted charm component (hatched 
histogram). Shown are the (a) decay length significance in the jet with the exclusive D*+ candidate; 
(b) decay length significance in the jet opposite the exclusive D*+ candidate; (c) distribution of the 
neural network based on jet-shape variables, and (d) distribution of the hemisphere charge. 



13 



hadronic decays of the Z° is found to be ng^cc = 0.0238 ±0.0048. Using Monte Carlo simulation 
the fraction of events from this source in the selected sample of events is estimated and, after 
normalising to the measured rate, subtracted from the tagged sample of D*+ events. The total 
contribution of these events to the exclusively tagged sample is found to be (0.2 ±0.1)%, where 
the error quoted is only due to Monte Carlo statistics. 

4.1.5 Systematic Errors of the Flavour Separation Method 

A number of errors are introduced by the flavour separation method employed. The errors 
quoted are relative errors on the charm fraction in the sample. 

• Detector resolution: The influence of the detector resolution on the tagging variables is 
studied in Monte Carlo by varying the resolutions in the tracking system by ±10% relative 
to values that optimally describe the data. The analysis is repeated and efficiencies are 
recalculated. The error quoted is the largest observed deviation in this study, and amounts 
to 0.6%. 

• Background modelling: The distributions in the tagging variables for background events 
are taken from estimator distributions determined using data. Possible differences be- 
tween the estimator and the true background distributions are studied using Monte Carlo 
simulations. The separation is repeated with the background distribution taken from 
Monte Carlo, and the differences are used as systematic errors. Similarly, errors in the 
determination of the normalisation of background in the sample will change the recon- 
structed b-fraction (and, since they are totally anti-correlated, the c fraction). This has 
been studied by varying the background within its total error. In total the error from 
these sources amounts to 1.7%. 

• Detector response: Possible inhomogeneities of the detector response as a function of cos 9 
are studied by repeating the flavour separation in bins of cos^, and by comparing the 
results with the overall determination. No signiflcant differences are found. 

• Hemisphere correlations: Part of the flavour separation is done in the hemisphere opposite 
to the reconstructed D*+ mesons. Small correlations are expected to exist between the 
two hemispheres, which possibly might bias the measurement of the flavour composition. 
In the flavour separation these biases are taken into account by calculating the tagging 
efficiencies in bins of xy)*+ as measured from the exclusive D candidate. The size of the bias 
is estimated by recalculating the tagging efficiencies in only one bin of 0.4 < xd*+ < 1.0, 
and repeating the analysis. This is done for all three tagging algorithms. The resulting 
error is 0.4%. 

• Charm modelling: The jet-shape analysis is sensitive to the modelling of the response to 
charm events, which is taken from Monte Carlo simulation. Possible modelling problems 
are investigated by comparing the network output distribution in an unbiased sample of 
hadronic Z° decays in data and Monte Carlo. All observed differences are assumed to 
come from charm modelling problems, and a systematic error of 1% is calculated. The 
same method was used in 0. 

• Charm and bottom multiplicity: The vertex finder employed is sensitive to the charged 
multiplicity from charm and bottom hadron decays in the sample. The multiplicity for 
heavy flavour decays in the Monte Carlo has been varied by reweighting simulated events. 
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corresponding to the current experimental bounds of ±0.2 tracks for charm decays, and 
±0.35 tracks for bottom decays Similarly the hemisphere charge technique is sensitive 



to the multiplicity, and its error is estimated using the same procedure. Overall this results 
in an error of 0.6% for charm and 0.5% for bottom. 

B hadron lifetime: The B hadron lifetime has been varied within its current experimental 
limits: In the D hemisphere the lifetimes of the different B species have been varied inde- 
pendently by ±0.07 ps for the B+, ±0.08 ps for the B°, and ±0.12 ps for the Bs [0. In the 
hemisphere opposite to the D meson, the mean B hadron lifetime of (1.549 ± 0.020)ps 
has been used and changed within its error. The total error is found to be 0.7%. 

Charm lifetime: The lifetimes of the weakly-decaying charmed hadrons D° and D"*" has 
been varied independently by ±0.004 ps for the D'^, and ±0.015 ps for the D+ ||22|| , 
corresponding to a total error of 0.4%. 

mixing: Mixing of neutral B mesons changes the correlation between the charge 
of the quarks in the two hemispheres of the event. The hemisphere charge method is 
sensitive to the assumed value of the mixing through the bottom tagging efficiency. Two 
different mixing parameters have to be considered: the average B mixing in the hemisphere 
opposite the tagged D*+, and the effective mixing in tagged D*"*" events. The latter is 
different because the mixture of B mesons in D*+ tagged events is different from an 
unbiased sample. The derivation of this effective mixing is described later in this paper 
in section |6.2.2| , and is determined with an error of ±17%. Compared to this the error 
on the average mixing is very small and is neglected. Applied to the b/c separation this 
translates into an error of 0.3%. 

Fragmentation model: The fragmentation functions used in Monte Carlo simulation for 



heavy flavour events influence all three tagging algorithms. The Peterson scheme 
is used for bottom and charm hadrons. To estimate the influence on the results the 
mean scaled energy of primary charm and bottom hadrons has been varied around their 



measured values of (xb) = 0.702 ±0.008 and (xd*+) = 0.510 ±0.009, as suggested in ||l6 
and the separation has been redone. This results in an error of 0.6%. 

The flnal charm fraction in the D*+ tagged sample for the selected xd*+ range is found to be 

/D*"" = 0.774 ± 0.008 ± 0.022 , (9) 
where the flrst error is statistical, the second one systematic. 



4.2 The Inclusive Charm Tag 

The second, more inclusive, charm tag relies on the very special kinematical properties of the 
decay D*"*" D°7r+. Because of the small mass difference of only 145 MeV between the D*+ 
and the D° very little phase space is left for the pion. In the laboratory frame this pion, called 
the "slow pion" in the following, is emitted essentially in the direction of the D*+ meson, with 
a maximal transverse momentum pt relative to the D*+ direction of flight of 39 MeV. A charm 
tag is constructed from this by looking for an enhancement in the density of tracks along the 
D*+ flight direction. 

The fraction of slow pions from cc events is enhanced by requiring: 
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• 1.0 GeV <p^< 3.0 GeV , 



where is the momentum of the pion candidate. Kaon and electron contamination in the slow 
pion candidate sample is reduced by using the particle identification power of the drift chamber 
in the OPAL detector, requiring 

• ^dE/dx > 0.02 , 

• n^E/dx > 20 , 

where W^dE/dx ^^e dE/dx probability for a pion, and UdE/dx the corresponding number of 
measurements used, as defined in section ^A\ . 

The fiight direction of the D*"*" meson is reconstructed inclusively by an iterative proce- 
dure. It uses the fact that decay products from heavy mesons are on average harder and more 
coUimated than those from fragmentation tracks, leading on average to higher values of the 
rapidity^ y = {l/2)\og[{E + pz)/{E — Pz)] with respect to fiight direction of the D*+ meson. 
The decay products are selected by first grouping all tracks and unassociated neutral clusters 
in the event into jets using the cone algorithm, described in section |^. The jet axis is com- 
puted from all particles in the jet after removing the slow pion candidate itself. If tracks or 
clusters exist which have a rapidity measured relative to the direction of the jet, of less than 
2.5, the one with the smallest rapidity is removed from the calculation, and the direction is 
recomputed. This procedure is repeated until all particles have a rapidity value above 2.5, or 
the number of tracks or clusters is less than two. In this case the original jet direction is used. 
The direction determined in this manner is used as an estimate of the D*+ fiight direction. The 
resolution, measured as the width of the pf distribution at 50% of its maximal value, is found 
to be cr{pt) = 0.056 GeV. The efficiency with which D*+ D°7r+ decays are selected using this 
method is around 40% in cc events, and around 20% in bb events. A similar procedure was 



first introduced in 25 



4.3 The Bottom Tag 

A pure sample of bottom events is selected using a lepton tag. This sample will be used in the 
measurement of f (b — D*+X) as the fiavour tagged sample, and takes the place of the exclusive 
D*+ tag described earlier. Electrons and muons are identified as described in section [1.1| . The 
sample is purified by requiring that electrons have a momentum larger than 2.0 GeV, and a 
transverse momentum relative to the jet direction larger than 1.1 GeV. Electrons are only 
reconstructed in the central region of the detector, if the polar angle is below | cos^| < 0.715. 
Muons have to have a momentum above 3.0 GeV, a transverse momentum larger than 1.2 GeV 
and I cos^^l < 0.9. According to the Monte Carlo simulation this sample has a bottom purity of 
(89.90 ± 0.14)%. Of the remaining events 33% are from semileptonic charm decays, and 66% 
are misidentified leptons. 

4.3.1 Systematic Errors of the Lepton Tag 

The purity of the lepton identification is taken from Monte Carlo simulation. The following 
systematic errors have been investigated: 

^The momentum is measured relative to the jet axis. Charged particles are assumed to be pious, neutral 
particles photons. 
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Detector resolution: The resolution of the tracking part of the detector is varied by 10%, 
resulting in an error on the purity of 0.1%. 

Bottom fraction in T? decays: The fraction of Z° — > bb events in the Monte Carlo is 
re-weighted to the one measured by OPAL: Tj^^/rhad = 0.2175 ± 0.0022 0]. The error on 
the measurement is used to calculate the corresponding systematic error of 1.0%. 

Heavy flavour fragmentation: The fragmentation parameters in the Monte Carlo have 
been varied to change the mean scaled energy of charm and bottom hadrons around their 
experimental values of 0.702 ± 0.008 and 0.510 ± 0.009 respectively. The error on the 
bottom purity is 0.8%. 

Decay modelling: The momentum distribution of the lepton produced in a bottom or 
a charm decay influences the tagging efficiency. Following the recommendations in [|l^ 
this was studied by reweighting the distribution in the Monte Carlo to different models. 
Models used are ACCM, ISGW and ISGW**. The largest observed deviation is used as 
a systematic error, resulting in an error of 1.4%. 

Semileptonic branching ratios: The semileptonic branching ratios B(b ^ £) and B(c — ^ £) 
have been measured at LEP. The spectra in B-decay are determined also at lower energy 



machines. The values recommended in |jT6| are used, and the Monte Carlo is re-weighted 
to these measured values. Systematic errors are derived from the errors on the branching 
ratios. The error on the purity is 0.3%. 

• Hadronic background: Around 6% of the sample of tagged leptons are hadrons, which 
were misidentified. In the error of the mistagging rate has been determined to be 
9.3% in the electron sample, and 9.0% in the muon sample. This translates into an error 
of 0.6% of the bottom purity. 

The total systematic error of the bottom purity in lepton tagged events is found to be 2.1%. 
Note that the knowledge of the lepton reconstruction efficiency is not required in this analysis, 
as discussed in section ^ 



5 Production of D*+ Mesons in TP — > cc and bb 
Decays 

In this part the total production rates Fcc/Fhad ■ f (c ^ D*+X) B(D*+ Xfix^) B(D° ^ K-7r+) 
and Fbb/Fhad ■ f (b D*+X) B(D*+ D°7r+) B(D0 K-7r+) are determined. The results for 
f (c — > D*"''X) is used later in the analysis to determine the relative partial width Fcc/Fhad- In 
addition the mean scaled energy of D*"*" mesons in T? cc events, (xd*+)c is measured, and 
the total multiplicity of charged D*+ mesons in Z° decays is given. 

The analysis is performed using the sample of fully reconstructed D*""" mesons in the 3-prong 
decay mode. To minimise the number of D*+ mesons not observed the XD*+-range for this part 
of the analysis has been extended to xd*+ =0.2. The reconstruction method is similar to the 
one described in the previous section for the exclusive charm tag. An important difference 
however is the treatment of the partially reconstructed D*+ decays. While previously they 
have been treated as signal for the purpose of tagging the primary event flavour, they are 
background for the determination of the total rate of D*"*" meson production in this particular 
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channel. Therefore the background subtraction is rediscussed in some detail, and a method of 
treating the contribution from such partially reconstructed decays is introduced. 

5.1 Background Subtraction 

Backgrounds for the purpose of the measurement of the production rate of D*"*" mesons are 
combinatorial background and partially reconstructed decays of D*+ mesons. The former is 
determined with essentially the same procedure as described above, except that the combina- 
torial background, after having been normalised, is not simply subtracted, but is fitted using a 
simple parametrisation. The contribution from partially reconstructed D*"*" decays is measured 
in the data with a special procedure. In this part of the analysis no requirement is made that 
only one candidate be found in the channel, unlike that for the sample of D*+ mesons used in 
the exclusive charm tag. This different treatment is possible since the 3-prong decay is very 
clean, and the number of partially reconstructed events is small. 

The combinatorial background in the sample is determined as before from the reflected pion 
estimator. The AM distributions obtained using this estimator are normalised to the candidate 
AM distribution for 0.16 GeV < AM < 0.2 GeV. It is then parametrised using an empirical 
functional form 

f(AM) = A (^(AM + B (AM?)f , (10) 

with A, B and C free parameters determined in the fit described below. The number of back- 
ground events is determined in 16 bins of x-d*+ between a;D*+ = 0.2 and 1.0, in the signal region 
0.142 GeV < AM < 0.149 GeV. The mass difference distribution for the 3-prong sample only, 
with the result of the background fit superimposed, is shown in figure ^(a) for the Xd*+ range 
between 0.2 and 1.0. 

From Monte Carlo simulation about 8% of the signal, after subtracting the combinatorial 
background, actually come from partially reconstructed D*+ mesons. They are mostly products 
of the following decays: D° ^ K'ti+ti^ (4.3%), ^ K-K+ (1.8%), 0° ^ 7r-7r+ (0.7%), D° ^ 
7r~7r"'"7r° (0.5%) and D° K^i^u (0.3%). The numbers in brackets indicate the predictions 
from the Monte Carlo simulation for the contribution to the full D*+ sample from each source. 

Instead of using the Monte Carlo predictions the total contribution is measured in data in 
a simultaneous fit to the Mq and the AM distributions of all candidate D*+ mesons. The Mq 
distribution is examined for candidates where the mass difference AM is inside the tight signal 
region of 0.142 GeV < AM < 0.149 GeV, and no Mq cut has been applied. Contributions from 
partially reconstructed D*+ mesons are in general characterised by a peak in the AM distri- 
bution at the position expected for correctly identified D*+ mesons, but no peak-like structure 
in the Mq distribution at the position expected for true D'' mesons. Therefore the difference 
between the number of reconstructed 0*+ mesons as derived from the AM distribution and and 
from the Mq distribution can be used as a measure of the fraction of partially reconstructed 
decays in the sample. A slight complication arises from the decay D° K~K+, which is ex- 
pected to peak just below the nominal D° mass. Monte Carlo simulation is used to account for 
this. 

The number of D*"*" candidates is extracted from the Mq mass spectrum using a fit. The 
combinatorial background is parametrised by an exponential function. The signal function has 
two contributions: the first describes the true 3-prong candidates, and is constructed from two 
gaussian functions, motivated from Monte Carlo simulation studies, with the second having 
three times the width and the same mean as the first. The second function parametrises 
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Figure 3: (a) Distribution of the mass difference AM = M^, — Mq reconstructed in the decay D*""" — >• 
D^TT'^, K~TT~^. Superimposed is the background distribution obtained from the background 

estimator discussed in the text, and the result of the fit to this background estimator. The arrows 
indicate the selected signal region, (b) Mq spectrum of D*"*" candidates, with Mq cut removed, and 
an additional AM cut applied. Shown are the data (points with error bars), the result of the fit as 
described in the text, and the two components from satellite and from fully reconstructed D*+ mesons, 
as obtained in the fit. The arrows indicate the selected signal region. 



the contribution from partially reconstructed decays in the signal region. It is dominated by 
satellite decays which cluster below the expected D° mass of 1.865 GeV. As expected from 
the kinematics of this decay the mass distribution is approximately gaussian, with a significant 
tail towards smaller masses. It is parametrised by a gaussian convoluted with an exponential 
function. The decay constant in the exponential is fixed relative to the width of the gaussian 
function to the value obtained in the Monte Carlo. The other decays contributing to the 
partially reconstructed sample are described by an additional exponential function, added to the 
parametrisation of the satellite decay. This essentially adds a tail to the satellite function, which 
extends into the nominal D" mass region. Monte Carlo is used to estimate the contribution 
coming from the D'' K~K"'" decay, which is not described by the tail. The shapes of the 
different fit functions have been tested in Monte Carlo simulated events, and are found to 
provide a good description of the mass spectra. Because of the complicated fit function, and 
because the fraction of partially reconstructed decays varies only slowly with xd*+, this fit is 
done in four equal-sized bins of Xd*+ between 0.2 and 1.0, instead of the 16 used in the AM 
fits. The number of D*+ mesons is obtained by integrating the signal function over the mass 
window 1.79 GeV < Mq < 1.94 GeV. The fraction of partially reconstructed D*+ events in this 
bin is then calculated from the difference between the number of signal events determined with 
this fit, and the sum of the number of signal candidates over the appropriate xd«+ bins found 
in the fit using the AM method. The Mq spectrum for all candidates, with the results of the 
fit superimposed, is shown in figure 0(b). 

Monte Carlo studies indicate that this method reliably reproduces the number of partially 
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reconstructed 0*+ mesons in the sample. The total contribution from all sources is predicted to 
amount to (7.9 ±0.5)%, while the fit in the Monte Carlo sample measures this to be (7.8 ±2.2)%. 
In the data the same procedure gives the contribution from partially reconstructed decays to 
be (8.1 ± 1.7)%, in excellent agreement. In total 8497 candidates are found, of which 3750 ± 24 
are background events, where the error given is the statistical error from the fit. 



5.2 Flavour Composition and Fragmentation Fits 

The main contributions to the sample of tagged events are the same as described in section |4.1.^j . 



The principal method for determining the flavour composition is the same as was described 
for the exclusive D*+ sample. The goal of this analysis is the determination of the absolute 
rate of D*""" production in charm and bottom decays. Therefore the observed number of D*+ 
mesons needs to be corrected for the reconstruction efficiency which can be done reliably only 
in the 3-prong sample. In addition knowing the efficiencies allows to constrain the shape of 
the fragmentation function to a particular shape. Previous studies 0, || have shown that 
the efficiency corrected fragmentation function in charm decays can be described well by the 



function of Peterson et al. 15 



The flavour separation is done in the xd*+ range 0.2 < a;D*+ < 1.0, subdivided into 16 bins 
of xd*+. The charm and the bottom fragmentation functions are constrained to the Peterson 
shape by convoluting the analytical fragmentation function with the effects of the hadronisation 
as predicted by the JETSET 7.4 model. A simultaneous flt is then done using the mean scaled 
energy and the information from the flavour separation procedure. The normalisations for both 
bottom and charm decays and the Peterson parameters Ec and are allowed to vary in the flt. 

The efficiency of the D*+ reconstruction is calculated in bins of Xd*+ separately for Z° — >■ bb 
and for Z° cc events in the Monte Carlo simulation. The efficiency is essentially constant as 
a function of Xd*+, with a small step at a;D*+ = 0.5 due to the change in the dE/dx and the 
cosr cuts. Typically it is (25.0 ±0.6)% for xd.+ < 0.5, and (30.0 ±0.5)% for xd.+ > 0.5, with 
the bottom and charm efficiencies being very similar. 

Some of the tagged D*+ events are expected to come from gluon splitting. The OPAL 
measurement for the gluon splitting rate is used 0,|2^, and the shape is taken from the Monte 
Carlo prediction as was done in 0]. The result of the separation is shown in figure ^ where 
the total efficiency-corrected yield, the charm component, the bottom component and the part 
from gluon splitting are shown. 

Integrating the fitted fragmentation functions over the full xd*+ range the product branching 
ratio in charm events is found to be 

Wrhad ■ f (c ^ D*+X) B(D*+ ^ D°7r+) B(D° ^ K-7r+) = (1.041 ± 0.020) x 10'^ , 

where only the statistical error has been quoted. The P^r degree of freedom is 1.28. 

Though not directly needed in this analysis the same fit returns information about the 
production of D*+ mesons in —>■ bb events, the total D*+ multiplicity, and the hardness of 
the D*+ fragmentation function in charm decays. In bottom events the product branching ratio 
is determined to be 

Tbh/rhad ■ f (b ^ D*+X) B(D*+ ^ D°7r+) B(D° ^ K-7r+) = (1.334 ± 0.049) x 10"^ . 

The errors quoted are only statistical. The correlation between the total production rate in 
Z° — > cc and in Z° — > bb events is found to be —23%. Adding the predicted gluon component, 
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D* 



Figure 4: Efficiency corrected yield of D*+ mesons as a function of tlie scaled energy Xj^.+ j 
l/-^haddA''D.+/dxj3.+ , for all candidates (filled points with error bars) reconstructed in the decay 
D*+ D''7r"'",D'' K~7r+, and the charm component after flavour separation (solid line). The open 
points are the bottom component, and the dashed line represents the result of the fit. Also shown 
is the predicted contribution from gluon splitting events (hatched area). The reconstruction is only 
done for X£).+ > 0.2, as indicated by the solid vertical line. 



and correcting for the branching ratios B(D*+ D°7r+) = 0.683 ±0.014 and B(D'^ K 7r+) = 
0.0383 ± 0.0012 p^, the total multiplicity of D*+ mesons in hadronic Z° decays is found to be 



nzo^B*+x = 0.1854 ± 0.0041 . (11) 

In addition the shape of the fragmentation function allows the determination of the average 
a:D*+ between and 1 in charm decays to be measured as 

(a;D.+ )e = 0.515 ±0.002 . (12) 

Again only the statistical error is quoted. Note that this measurement does not include the 
effects from D* mesons produced in gluon splitting. 

5.3 Systematic Errors of the Measurement of D*+ Production 

A number of systematic errors are investigated in connection with the D*""" rate measurements 
in charm and in bottom events. Note that the errors are totally anticorrelated, and only the 
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ones for the rate in Z° — ^ cc events are given. Errors of the determination of the mean scaled 
energy of D*+ mesons are discussed separately at the end of this section. The first group of 
errors are due to detector effects, resolutions and Monte Carlo modelling: 

• Track quality cuts: The effects of the track quality cuts are investigated by comparing 
the efficiency for each selection cut in data and Monte Carlo. An error of ±0.6% has been 
found to be sufficient to cover observed differences. 

• Fraction of silicon hits: The resolution of secondary vertices depends on the fraction of 
tracks which use measurements from the silicon micro-vertex detector. The fraction in 
Monte Carlo events has been re-weighted to the one measured in data. An error of ±0.4% 
is derived from the statistical precision of this procedure. 

• dE/dx modelling: The calibration of the specific energy loss, dE/dx, has been compared 
in data and Monte Carlo in samples of identified particles. Samples of kaons and pions 
are selected in decays oi (j) —>■ K~K"'" and K*^ vr+vr" mesons without using dE/dx 
requirements, and the calibration of dE/dx is measured in the data. In addition D*+ 
mesons are reconstructed without dE/dx requirements, and the results are compared 
with those quoted in table ^ For the applied cut of 2% on the kaon weight an error 
of 1.1% of the total rate was found in those candidates where a cut was applied. This 
error contains a contribution from the measurement of dE/dx itself, mainly due to the 
calibration, and from the requirement of at least 20 hits for the dE/dx measurement. 
Since a dE/dx cut is only used for xd*+ < 0.5 this translates into a reduced error on the 
total rate in charm events of 0.8%, and in bottom events of 1.0%. 

• Mass resolution: The invariant mass resolutions in data and Monte Carlo for the decay 
selected have been compared. The Mq resolution in the Monte Carlo is 27.5 MeV, the 
one in data 27.9 MeV. Depending on the Monte Carlo sample used variations in the mass 
resolution from sample to sample of up to a few MeV are observed. These differences 
correspond to changes in the momentum resolution of the detector of approximately 10%, 
and translate into a systematic error on the efficiency of 1.0%. 

• Rate of partially reconstructed D*+: The total contribution from partially reconstructed 
D*+ mesons is measured in the fit to be (8.1 ± 1.7)% in the selected sample. Of these 
(4.6 ±0.8)% are reconstructed as coming from the satellite decay mode. The Monte Carlo 
simulation predicts the fraction of partially reconstructed D*+ decays to be (7.9 ± 0.5)%, 
of which (4.8 ± 0.4)% are from the satellite decay mode. The rest, (3.3 ± 0.3)%, are from 
a number of other decays, as discussed in section |0| . In general very good agreement 
is observed between the predicted and the measured fractions, both within the Monte 
Carlo simulation, and between data and Monte Carlo. The uncertainty of the method 
is estimated from the statistical precision of the fit, and from the error of the individual 
branching ratios contributing to the partially reconstructed signal in the Monte Carlo. In 
addition a contribution of 0.5% is included to account for the finite Monte Carlo statistics 
available for this study. In total the relative error of the rate of partially reconstructed 
D*+ mesons contributing to the signal is estimated to be 22% of the rate of partially 
reconstructed mesons, which contributes an error of 1.5% to the total rate measurement. 

• Background subtraction: The combinatorial background in the sample is subtracted based 
on estimators derived from data. The quality of the procedure has been studied in the 
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Monte Carlo simulation. Within the available statistics no significant deviations are 
found. The total difference is less than 1%, and an error of 1% is assigned to this source. 



g — ^ cc: The expected contribution from gluon splitting has been subtracted from the sam- 
ple, for xd*+ > 0.2. The mean value used is the one measured in [0,^ of 0.0238 ± 0.0048. 



The contribution from this process has been varied within this error, which results in an 
error on the total rate of ±1%. Monte Carlo studies indicate that the shape of the g ^ cc 
component is not very dependent on the particular Monte Carlo model used. Comparing 
JETSET and the Ariadne Monte Carlo model an error of ±1% is assigned to this 



source, resulting in a total error from gluon splitting of ±1.4%. 

Heavy flavour fragmentation: The total rate of D*+ mesons produced in Z° —>■ cc and in 
— > bb decays has been obtained by integrating the fitted fragmentation functions from 
xy)*+ = to 1. This procedure is subject to a number of systematic uncertainties: 

— fragmentation function: The Peterson fragmentation function has been used for the 
results quoted. The different fragmentation models of Collins and Spiller ET^ and 



Kartvelishvili p8[ have been used to estimate its influence. The largest difference 
in the fitted rate has been used as the systematic error from this source. The error 
found is 1% of the rate in charm events, and 4% in bottom events. 
As a cross-check the fits have been repeated using the QCD inspired fragmentation 



function from Nason et al. p9[. The results are compatible with the ones obtained 
using the fragmentation models. The error on the fitted parameters however also 
reflect the limited statistics available for the initial tune of the model using low 
energy data. It is therefore used as a cross-check rather than to give an additional 
error. 

— b decay modelling: A comparatively large fraction of b — D*+ decays are not ob- 
served since only events with x > 0.2 are considered in this analysis. In addition to 
the error from different fragmentation models discussed above, an error introduced 
by this extrapolation on the measured b-rate has been studied as in p by consid- 
ering the different types of decays contributing to the spectrum, and investigating 
the differences in the extrapolation introduced by each of the components. The 
differences found amount to an error of 2% on the b — > D* rate measurement. 

— excited D meson production: The shape of the fragmentation functions is influenced 
by the presence of D excited states in the decay chains. About 32% of all D*+ 
mesons have been measured as originating in decays of excited charm mesons. The 
contribution from these has been varied by ±18% around the mean value of 32%, as 
was done in p|. The resulting error is 0.2% on the rate. 

The total error from all fragmentation modelling issues is 3% for charm, 5% for bottom. 

In addition the errors discussed in the section on the flavour separation apply to the rate 
measurements as well. A complete breakdown of errors for the measurements is given in table 13. 



Only a few of the errors listed for the rate measurement have a significant effect on the 
determination of the mean scaled energy of D*+ mesons in charm decays. The most important 
error comes from the extrapolation of the fragmentation function into the unmeasured region 
below xj)*+ = 0.2. Making the same comparisons between different fragmentation models as 
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chaiiii 


bottom 






error source 






detector resolution 


track quality cuts 


0.006 


0.008 


0.0008 






fraction of silicon hits 


0.003 


0.004 


0.0004 




> 0.001 


dE/dx modelling 


0.008 


0.013 


0.0012 






total detector resolution 


0.011 


0.016 


0.0015 


0.001 


D*+ reconstruction 


mass resolution 


0.010 


0.013 


0.0013 


■ 




subtracting partially reconstructed decays 


0.016 


0.020 


0.0019 




> 0.002 


D*+ background subtraction 


0.010 


0.013 


0.0013 






g^ cc 


-0.016 


-0.019 


-0.0018 




0.004 


total D*+ reconstruction 


0.026 


0.030 


0.0032 


0.005 


flavour separation 


/b statistical error 


0.010 


0.013 








background modelling 


0.017 


0.023 








hemisphere correlation 


0.004 


0.005 








charm modelling 


0.010 


0.013 








bottom multiplicity 


-0.005 


+0.007 






> 0.002 


charm multiplicity 


+0.006 


-0.008 








bottom lifetime 


-0.007 


+0.009 








charm lifetime 


+0.003 


-0.004 








B mixing 


0.003 


0.004 




y 




fragmentation model 


0.012 


0.053 


0.0041 




0.007 


b decay model 




0.027 


0.0021 






excited D meson production 


0.002 


0.003 


0.0003 




0.002 


total flavour separation 


0.030 


0.070 


0.0047 


0.008 


total 


0.040 


0.078 


0.0059 


0.009 



Table 3: List of systematic errors relevant in the determination of the rate of D*+ mesons produced 
in TP cc and in TP bb events, the total multiplicity of D*"*" mesons in TP decays, and the mean 
scaled energy. For the error of the mean scaled energy only errors which are larger than 0.001 x 10~^ 
are listed separately, otherwise the total contribution from a class of errors is given. A sign in front 
of an error indicates the direction of change under a positive change of the variable. The last three 
errors from the flavour separation include the contribution from the same sources through the D*+ 
reconstruction. Note that all flavour separation errors are anti-correlated between the charm and the 
bottom result. 
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described above a modelling error of ±0.007 for the mean scaled energy has been determined. 
The uncertainty in the modelling of gluon splitting results in a further error of 0.004. Effects 
related to the flavour separation have much less of an effect on the mean scaled energy. Taken 
together all other errors except the above-mentioned modelling issues contribute another 0.004 
to the error. The reconstruction of D*""" mesons, in particular the background subtraction, 
contribute an additional overall error of 0.002. 

The final results, including all systematic errors, for the hadronisation fractions of charm 
and bottom quarks into D*""" mesons are found to be 

Tcc/rhad ■ f (c D*+X) B(D*+ ^ D°7r+) B(D° ^ K-7r+) = (1.041 ± 0.020 ± 0.040) x 10^^ , 

and 

Tbh/rhad ■ f (b ^ D*+X) B(D*+ ^ D°7r+) B(D° ^ K-7r+) = (1.334 ± 0.049 ± 0.078) x 10"=^ , 

where the first error is statistical, the second systematic. The total multiplicity of D*+ mesons 
in hadronic Z° decays is found to be 

nzo^B*+x = 0.1854 ± 0.0041 ± 0.0059 ± 0.0069 . 

The last error quoted is due to external branching ratios. From the shape of the fragmentation 
function the average mean scaled energy xd*+ of D*+ mesons in charm decays is determined to 
be 

(xd.+ )c = 0.515 ± 0.002 ± 0.009 . 

The errors quoted are statistical and systematic, respectively. The results quoted is for the 
primary production of D*+ mesons in Z° — > cc events, and does not contain contributions from 
gluon splitting events. The results are in agreement with other measurements at LEP and with 
the previous OPAL determination |l7|,^[3l[]. 



6 Measurement of f (c ^ D*+X), f (b ^ D*+X) 
and of Tcc/rhad 

The hadronisation fraction f (c — » D*+X) is measured by the simultaneous detection of charm 
quarks in two jets of the event. In one jet, a charm meson candidate is identified by using the 
exclusive D*"^ meson tag discussed in section [4.1| . In the other jet the inclusive tag presented in 
section is used. Background is suppressed by requiring that both D*+ tags are of opposite 
charge where the charge of the D*"*" is given through the charge of the slow pion candidate 
from the D*+ D^n^ decay. By determining both the number of exclusively reconstructed D 
mesons and of simultaneously tagged jets, the hadronisation fraction can be calculated. 

The hadronisation fraction f (b ^ D*+X) is measured in an analogous fashion by the simul- 
taneous detection of a lepton in one jet and an inclusively reconstructed D*+ meson in another 
jet of the event. The exclusive charm tag is replaced by a lepton tag, optimised for the selection 
of bottom events. The inclusive tag based on the slow pion is used in the opposite jet. 

Combining the hadronisation fraction f (c ^ D*+X) with the measurement of the production 
of D*+ mesons in Z" —>■ cc events the relative charm partial width Pcc/rhad is calculated, as 
discussed at the end of this section. 
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6.1 The Tagged Samples 



Two samples are used in this analysis, one tagged by the presence of an exclusively reconstructed 
D*+ meson, the other tagged by a lepton. Both are selected as described in section ^ After 
all cuts a sample of 27005 D*+ candidates has been found, of which 11 366 ± 107 are estimated 
background events. Of the selected D*+ events a fraction of 0.774 ± 0.023 are reconstructed 
as coming from TP — ^ cc events. The number of tagged electron and muon candidates is 
determined to be 43 579, of which 4445 ± 64 do not originate in bottom decays. 

Slow pions are sought in the secondary jet in the samples of events tagged by either a 
D*+ meson or a lepton. The reconstruction of slow pion candidates proceeds as described in 
section W?l. Event samples are prepared for the correct charge combination D*"''-7r~ and for any 



charge combination £-7r, respectively. Background in the slow pion sample is estimated from a 
D*"-7r~ sample, prepared using a selection in a sideband of AM > 0.18. A clear enhancement 
is visible at low values of pi as shown in figure ^(a) for the D*+ tag, and in figure ^(b) for the 
lepton tag. A much smaller enhancement is visible in the sideband selected double tag, shown 
in figure |](c). The number of double tagged candidates, counted below pi < 0.01 GeV^, is 2146 
D*+-7r"'" candidates, and 5502 i-ir candidates. 

The efficiency to reconstruct a slow pion in the presence of a fully reconstructed D*+ meson 
or a lepton in the other jet is determined in Monte Carlo simulated events. The slow pion is 
sought in the secondary jet of the event, and its reconstruction efficiency is calculated. This 
procedure is done individually in each of the five exclusive D*+ modes, and for the lepton tagged 
sample. The final efficiency is calculated by reweighting the Monte Carlo efficiencies to reflect 
the mixture of tagged events in the data. After all cuts the efficiencies are found to be 

e^,^ = 0.384 ±0.004 and ej^*^ = 0.189 ± 0.006 , (13) 

for c —>■ D*+ and b —>■ D*+ events, respectively. In lepton tagged events the efficiencies are 

e^^ = 0.414 ± 0.029 and e^^ = 0.193 ± 0.005 . (14) 

The errors quoted are purely statistical. 



6.2 Composition of the Double Tagged Sample 

A number of different classes of events contribute to the sample of double tagged events. At very 
low p1 a significant fraction of the candidates are due to slow pions from the decay D*+ D^vr^, 
both in charm and in bottom decays. The signal in bottom decays, while similar to the one 
in charm, has a broader distribution in pi. Background in the sample comes from a number 
of different processes. The dominant source is random tracks that pass the applied cuts. This 
combinatorial background falls significantly more slowly with increasing pi than does the signal, 
and does not exhibit the characteristic enhancement at very low pi- 

A small but important background is slow pions from fake double tag candidates. A double 
tagged event is denoted a "fake double tag" , if the slow pion candidate found in the inclusive 
tag is correctly identified, but the fully reconstructed D*+ meson or the lepton in the other jet 
is wrongly identified. Such events contribute to the peak in the pi spectrum, and need to be 
subtracted from the sample. 

In the following each of the different parts of the candidate distribution will be briefly 
discussed. In the last part of this section the method used to count the number of double tags 
from charm decays is presented. 
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6.2.1 Combinatorial Background and Fake Estimation 

The combinatorial background is the dominant background source. Its shape is estimated using 
events with the wrong charge correlation between the D*+ and the slow pion in the opposite 
hemisphere, which in addition are reconstructed in a sideband of AM on the exclusive side, be- 
tween 0.18 (0.19) GeV < AM < 0.20 (0.22) GeV (numbers in brackets are for the semileptonic 
channels). Searching for slow pions in the secondary jet relative to the exclusive candidates 
tagged in the sidebands, a signal for charm production is observed at low (figure ^(c)). This 
signal has two contributions, one from fake double tags, another small one from incompletely 
reconstructed D*"*" meson decays in the D-jet. In both cases a true slow pion is found in the 
secondary jet. The total fraction of fake double tags in the background is measured in the 
sidebands, and is subtracted from the total number of double tags, as described in |673| . 



6.2.2 Contribution from Bottom Events 

The shape of the pf signal in bottom events is determined in data from the lepton-slow pion 
double tagged sample, which is about 90% pure in bottom decays. The fraction of events in the 
D*+-7r^ double tagged sample originating from bottom decays is determined from the known 
fraction of b-events in the D*+ sample, and the efficiency to tag a slow pion in a b-decay in the 
secondary jet. 

The situation is slightly complicated by mixing in the neutral B system. If mixing has 
occurred in either hemisphere, the charge correlation between the primary quark and the cor- 
responding D*+ mesons is changed, and the correlation between the charge of the slow pion 
track and the fully reconstructed D*+ candidate is opposite to the unmixed case. This fraction 
of events migrates out of the signal sample in this analysis. The total probability in bottom 
events that mixing has destroyed the charge correlation is given by 

XeS^ = XD*+(l-X7r+„ow)+X7r+„ow(l-XD-+) (15) 

where Xtt+sIow' Xd*+ are the effective mixing parameters applicable to the slow pion and the D*+ 
sample, respectively. The two mixing parameters are equal, since a D*+ tag is used in both jets. 
The majority of D*+ mesons in Z° — > bb events originate from decays of the B° meson. This 
fraction is estimated using semileptonic B decays to be = 0.790lo;i2 [Hi- ^ small percentage 
of D*"*" mesons are also expected from Bg mesons, which also mix. The number of D*+ from Bg 
mesons was estimated to be Vs = 0.033 ± 0.015. The average mixing in the neutral B system 



is determined from the world average value for the mixing parameter, Xd = 0.175 ± 0.016 P2 



For the mixing parameter of the Bg meson the current world average limit of Xs > 0.49 at 95% 
confidence level has been used. In this analysis Xs is varied between 0.49 and 0.50. The 
effective mixing seen by the D*+ mesons is given by 

Xd*+ =rdXd + Ts Xs- (16) 

In addition D*+ mesons with the wrong sign can be produced in bottom decays, when a c quark 
is produced in the decay of the W. This can be expressed in terms of a mixing-like parameter 
(y)- As in [0 a value of Cd = 0.025 ± 0.025 is used. The effective mixing parameter for the 
D*+-7r^ double tag sample is estimated to be 

X^s^ = 0.289 ± 0.050. (17) 

This number is in agreement with a direct measurement of the effective D*+ mixing in p!8|. 
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6.3 Determination of the Number of Double Tags 

The number of double tagged events in the sample is estimated from the three p"^ distributions 
by a simultaneous fit. The right sign sample is fitted as a superposition of true signal from 
charm and bottom decays, a contribution from fake double tagged events, and background. 
Each of the individual contributions is described by an exponential function 

F{pf)j=bgd,h,c = aj exp{bjPt) . (18) 

The parameters of the signal originating from bottom decays is determined in the i-n double 
tagged sample. The fake distribution is measured in a double tagged sideband sample, where 
the exclusive candidates are selected in a sideband of AM. Combinatorial background is fitted 
for in the different distributions. 

All three distributions are simultaneously fitted using the following parametrisations. The 
shape of the lepton tagged spectrum is parametrised by 

= F(p?)bgd + fl f (b - D*+X) el + f (c ^ D^X) , (19) 

where the purities and have been given in section and e^, e}^ are the efficiencies to 
find a slow pion in the secondary jet in the presence of a lepton in the other jet, in bb and 



cc events respectively, as quoted in section |6A| . Since both relative signs between leptons and 



pions are used eq. |T9] does not depend on the mixing in the lepton tagged sample. The right 
sign signal distribution is described by 

i^(P?)D.. = i^(p')bgd + F(p?)f^,^ 

+ (1 - xSr") (1 - /c°") f (b - D-X) F{pI)^ el,^ (20) 
+ /r"f(c-D*+X)F(p?)ee^.^. 

The charm fraction /°*^ fulfils the condition /°*^ = (1 — /b**^), and ej^* ^, Cd* tt the 
efficiencies to find a slow pion in the presence of a D*+ meson in the other jet in cc and bb 
events respectively. The mixing probability x^r^ has been determined in section |6.2.2 . 



The contribution from fake double tags in the D*+ — tt^ double tagged sample is measured 
in the sideband tagged sample, as described above. The pi spectrum in the sideband tagged 
sample is parametrised by 

i^(p?)side = « (i^(p?)bgd + i^(p?)fake) , (21) 

where -F(p^)fake contains contributions from fakes in the double tagged sample and from partially 
reconstructed D*""" mesons. The contribution from fake double tags is assumed to have the 



fiavour composition as the real signal, as given by the two last lines of eq. and the same 
functions are used for -F(p^)fake as for signal events. The absolute contribution from the fakes is 
obtained by rescaling the fitted fake contribution by the ratio a of the number of background 
candidates in the sideband sample to that in the signal sample. 



Equations |T9|, 20 and |2T| are fitted simultaneously. Free parameters in the fit are the hadro- 
nisation fractions f (c — D*+X) and f (b — > D*+X), the normalisation of the fake rate, the 
shape parameters a and h of bottom and charm slow pion signals, and the background pa- 
rameters. The different spectra are illustrated in figures |^(a) to (c), with the results of the fit 
superimposed in each case. The number of double tagged events in the D*"''-7r~ double tagged 
sample in Z° cc events, A^p,+ ^_, and in the £-7r tagged sample, iV^"^, are determined by 
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Figure 5: Spectrum of the squared transverse momentum of double tagged candidates after applying 
all cuts of the (a) signal candidates, reconstructed with opposite charges in both jets. The non-charm 
component is indicated by the hatched area; (b) lepton-tagged candidates. The points are the data 
and the lines the results of the fit; (c) background candidates, reconstructed in AM sidebands. 
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integrating the fitted signal functions in tlie double tagged samples between = GeV^ and 
p'f = 0.01 GeV^. They are found to be iV^*^ = 702 ±44, over a background of 1444 ± 18 events, 
and Nf^ = 934 ± 80 over a background of 4568 ± 45 events, respectively. Correcting for the 
branching ratio 3(0*+ —>■ D^tt^) = 0.683 ± 0.014 the hadronisation fractions are found to 
be 

f(c^D*+X) = 0.222 ±0.014 
f(b^D*+X) = 0.173 ±0.016, 

where the error quoted is only statistical. 

Combining the measurement of the hadronisation fraction f (c — D*'''X) with the result for 
the total rate of D*+ mesons in Z° cc events the charm partial width relative to the total 
hadronic width of the Z" is found to be 

Tcc/rhad = 0.180 ± 0.011 , 
where the error quoted is purely statistical. 

6.4 Systematic Errors 

In this section the different sources of systematic errors for this part of the analysis are discussed. 
A full breakdown of all errors considered is given in table ^. 

• Detector resolution effects: The effects of detector resolution modelling have been dis- 
cussed in section |4.1.5| . The total error applicable from these sources to this measurement 
is ±1.8%, with the dE/dx error being dominant. 

• Flavour separation: The systematic errors of the flavour separation determined in sec- 
tion |17T]3| are used to determine the corresponding systematic errors on the hadronisation 



fractions. Many of the errors for Fcc/Fhad are correlated to that of the charm hadronisa- 
tion fraction. This correlation is taken into account in calculating the final errors. 

• Lepton identification: The systematic errors discussed in connection with the lepton 
identification are applicable to the measurement of f (b D*'*'X). 

• D*+ reconstruction: Fcc/Fhad depends on the total rate of D*+ in events. All errors 
discussed in section | are applicable to this measurement as well. 

A number of additional errors are introduced through the inclusive charm tag: 

• Heavy flavour fragmentation: The efficiency for flnding slow pions in Z° cc events 
has been calculated in Monte Carlo. Systematic errors from the modelling of the heavy 
flavour fragmentation are studied as described in section |5.3| . The total error from this 
source is ±1.4% on the slow pion efficiency. In table ^ this error is combined with the 
fragmentation error from the b/c separation. 

• Excited D meson modelling: The influence of excited D meson production on the slow 



pion efficiency is studied as described in section p.3| , resulting in an efficiency error of 
0.1%. 
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Source 


f (b ^ D*+X) 


f (c ^ D*+X) 




detector resolution 


track quality cuts 


0.0000 


0.0000 


0.0011 


dE l^x modelling 


0.0017 


0.0025 


0.0018 


detector resolution 


0.0010 


0.0013 


0.0009 


flavour separation 


/b stat error 


0.0002 


0.0022 


0.0025 


background modelling 


0.0003 


0.0030 


0.0028 


hemisphere correlation 


< 0.0001 


0.0007 


0.0007 


charm modelling 


0.0002 


0.0018 


0.0017 


bottom multiplicity 


< 0.0001 


-0.0009 


-0.0008 


charm multiplicity 


0.0001 


+0.0011 


+0.0010 


bottom lifetime 


0.0001 


-0.0012 


-0.0011 


charm lifetime 


0.0002 


+0.0011 


+0.0008 


fragiiieiitatioii riKKloUiiig 


< 0.0001 


0.0011 


0.0018 


lepton reconstruction 


lepton fragmentation 


0.0030 


< 0.0001 




bottom fraction in Z° decays 


0.0017 


< 0.0001 




lepton decay model 


0.0050 


< 0.0001 




semi-leptonic BR 


0.0010 


< 0.0001 




le])t()ii fake rate 


0.0020 


< 0.0001 




D*+ reconstruction 


mass resolution 






0.0018 


subtracting partially reconstructed decays 






0.0027 


D*+ background subtraction 






0.0018 


double tag, slow pion reconstruction 


slow pion efficiency 


0.0040 


0.0010 


0.0025 


excited D meson production 


0.0040 


0.0010 


0.0008 


fitting procedure 


0.0081 


0.0104 


0.0085 


gOgo niixing 


0.0002 


0.0016 


0.0011 


jet resolution 


0.0024 


0.0031 


0.0025 


jet-jet correlation 


0.0012 


0.0015 


0.0012 


g ee 


< 0.0001 


-0.0003 


-0.0021 


total error 


0.0121 


0.0133 


0.0123 


external inputs 


B(D*+ ^ D°7r+) 


0.0034 


0.0044 




B(D0 ^ K-7r+) 






0.0063 



Table 4: List of systematic errors contributing to f (c -^^ D*+X), f (b -^^ D*"''X) and Tcc/rhad- A sign 
in front of an error indicates the direction of change under a positive change of the variable. A detailed 
explanation of the different errors can be found in the text. 
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• Jet definition: The resolution of tlie signal is dominated by the resolution of the 
direction of the jet axis. Possible modelling problems in the Monte Carlo have been 
checked by comparing the number of charged tracks and of neutral clusters used in the 
calculation of the jet axis. The Monte Carlo distributions have been re-weighted to the 
data distributions, and the efficiencies are re-determined. An error of 0.9% in charm, 
and 3% in bottom events has been found. As a cross-check the fit has been repeated 
fixing the fitted resolution^ of the distribution of 0.056 GeV to its Monte Carlo value 
of 0.061 GeV. Consistent results have been found within the quoted errors. 

• Jet-jet correlation: The efficiency for identifying a slow pion is influenced by the presence 
of an exclusively reconstructed D*+ meson in the other hemisphere. This bias is taken 
into account by calculating the efficiency in Monte Carlo in events where a D*"*" meson 
has been reconstructed in the other hemisphere. A small error remains if the energy 
distribution of the secondary jet is not modelled properly in the Monte Carlo. This has 
been evaluated by reweighting the Monte Carlo energy distribution to that observed in 
the data, and recalculating the efficiency. The resulting error is 0.7% of the efficiency. 

• Background subtraction and fitting procedure: Systematic effects introduced by the fitting 
procedure for background determination in the double tagged sample have been studied 
in Monte Carlo simulation. 

— The complete fit has been repeated in Monte Carlo. The number of double tags 
reconstructed is well reproduced within its statistical errors. The difference is 0.5%, 
and this is used an additional systematic error. 

— The modelling of the combinatorial background in the fit has been tested in the 
Monte Carlo by repeating the fit with the true Monte Carlo background. The fit 
has been done by either constraining the shape of the background to be the same 
in all double tagged samples, or by fitting it individually in each one. The observed 
difference of 2.5% has been assigned to this source. 

— The shape of the bottom component is estimated from the lepton sample. Possible 
biases have been investigated in the Monte Carlo by repeating the fit with the true 
bottom signal. The results agree to within 2%, and this has been used as the error. 

— The number of fake double tags predicted by the fit in the Monte Carlo has been 
compared to the known number of fake double tags. They agree to 10%. The number 
of fake double tags in the data has been measured in the fit to be (8.5 ± 2.2)% of 
the combinatorial background. Monte Carlo predicts this to be (9.9 ± 0.8)%, which 
is compatible within errors with the measured rate. The error from this source 
is estimated by varying the fitted fake rate within its error, and by additionally 
assigning the difference between data and Monte Carlo as an systematic error. In 
total this gives an error of 3.5% of the signal in charm events, 3.0% in bottom events. 

— B-mixing: The uncertainty due to mixing in the neutral B sector has been studied 
by varying the effective mixing parameter Xes within its error, or 0.5% of the final 
result. 

— Fit function: As a cross check the analysis has been repeated with different parametri- 
sations for both the signal and the background distributions. The background 

^The resolution is defined as the width of the distribution at 50% of its maximum value. 
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is modelled using a polynomial {a + b + c p^)"^, and no differences are found 
for the results. The signal function has been replaced by a gaussian-like function 
exp(— (a + b pfY), which also results in consistent results. 

• g — > cc: The contribution to the tagged samples from gluon splitting events has been 
discussed in section ^.3| . It has been re-evaluated for the inclusively tagged jet, where 
Monte Carlo simulation predicts this to be (0.64 ± 0.08)%. A systematic error of 0.12% 
has been assigned to this source. 

A list of all systematic errors is given in table ^ 



7 Results and Conclusions 

A double tagging technique has been used to measure the hadronisation fractions of charm and 
bottom quarks into charged 0*=*= mesons. They are determined to be 

f(c^D*+X) = 0.222 ± 0.014 ± 0.014 ± 0.004 , 
f(b^D*+X) = 0.173 ± 0.016 ± 0.012 ± 0.003 , 

where the first error quoted is statistical, the second systematic, and the third one due to 
external branching ratios. 

From the number of D** mesons observed in the decay D*+ — D*^7r+,D*^ K^vr^ the 
multiplicity of D*"*" mesons in hadronic Z decays is measured to be 

nzo-.D-+x = 0.1854 ± 0.0041 ± 0.0059 ± 0.0069 . 

Applying bottom tags based on lifetime, jet shape and hemisphere charge information in the 
event the charm and the bottom components have been separated, and the individual produc- 
tion rates are found to be 

Tcc/rhad f (c ^ D*+X) B(D*+ ^ K"7r+7r+) = (1.041 ± 0.020 ± 0.040) x 10"^ , 
Tbh/rhad f (b ^ D*+X) B(D*+ ^ K-7r+7r+) = (1.334 ± 0.049 ± 0.078) x 10"^ , 

The mean scaled energy of D*+ mesons in charm events is determined from the fragmentation 
function to be 

(xd*+)c = 0.515 ± 0.002 ± 0.009 . 

From the hadronisation fraction f (c ^ D*+X) and the total production rate of D*+ mesons 
in Z° — > cc events, Fcc/Fhad is found to be 

Tcc/rhad = 0.180 ± 0.011 ± 0.012 ± 0.006 . 

Here the first error is statistical, the second one describes internal systematics, and the last one 
is due to branching ratio B(D° — K~n~^) = 0.0383 ± 0.0012 ||2^. The correlations between the 



rate measurement and the hadronisation fraction are taken into account in this calculation. A 
detailed breakdown of the systematic error is given in table ^. The measurements presented 
in this paper are based on the full data sample of almost 4.4 million events collected with the 
OPAL detector at LEP at a centre-of-mass energy of about 91 GeV. For the first time the charm 
partial width has been measured without significant input from lower energy experiments, and 
in particular without assumptions about the centre-of-mass energy dependence of heavy flavour 
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fragmentation. In a previous OPAL publication the hadronisation fraction f (c — D*+X) was 
derived from measurements at lower energy e+e" machines to be f (c D*+X)iow energy = 
0.262 ± 0.019 ± 0.010, where the last error is from the branching ratios B(D*+ D^tt"*") and 
B(D° — > K^TT^). Assuming that the only correlation between the low energy result and the 
result presented in this paper is from the branching ratio B(D* D^vr^), the low energy 
number is 1.4 standard deviations higher than the OPAL one, which is compatible with the 
assumption that the sources of D*""" mesons at lower energies and at LEP energies are the same. 
If the low energy hadronisation fraction were to be used instead of the one measured by OPAL 
Fcc/Phad would be lower by 15%. 

In conclusion good agreement is found with the prediction of the Standard Model of 
Pcc/Phad = 0.172, and with other measurements of the same quantities at LEP 
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